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PT distributions) are sensitive to the geometry (the angles between jets) and the 
colour structure of the underlying hard process (topology of the three-prong parton 
antenna). The results and techniques presented here could be extended to other 
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1. Introduction 

Multi-jet ensembles are typical for hard hadron-hadron coUisions. Production of 
large-pt hadrons has mainly the topology of a four-jet event; in W^, or large-pt 
photon production there arc three hard jets: two formed by the initial state radiation 
and one from the hard parton recoiling against the electro-weak boson. 

As a first step in improving the QCD description of processes involving more 
than two jets in the final state, we study three-jet hadronic systems in a simpler 
environment — without hadrons in the initial state. In this paper we continue the 
QCD analysis [1] of near-to-planar three-jet events in e'^e~ annihilation into hadrons 
in the regime 



where T,Tm, and are the thrust, the thrust major and the thrust minor, respec- 
tively. 

We study the distribution in Kout in three-jet events and aim to reach the same 
theoretical accuracy that is available these days for the description of the typical e'^e~ 
hadronic systems, that is two-jet events. The coUinear and infrared safe (CIS) dis- 
tributions in thrust T, C-parameter, jet masses and broadenings B in the region 
of two narrow jets (1— T, C, M'^/Q'^, 5^1) have been intensively studied in recent 
years. Both the perturbative (PT) contribution and the leading non-perturbative 
(NP) corrections have been computed at high accuracy. The PT results [2] involve 
all-order resummation of double- (DL) and single-logarithmic (SL) contributions and 
matching of the approximate resummed expressions with the exact second order ma- 
trix elements. The leading NP corrections [3]- [8] have been computed at two- loop 
order to take into account effects of non-inclusiveness of jet observables. The 1/Q- 
suppressed NP power corrections arc needed not only to make quantitative predic- 
tions of event shapes but to expand our knowledge of hadronization phenomena and 
of the interaction in the confinement region in general. 

In [1] we have computed the SL-rcsummed PT expression for the integrated dis- 
tribution S(/Cout) (the ii'out-spectrum integrated up to a given i^out) in the three-jet 
region (1.1). We considered the cumulative out-of-plane momentum both in the total 
phase space (-f^^t) restricted right region (K^^), that is in the narrow- 

jet hemisphere. The results turned out to be quite interesting and informative. In 
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particular at the SL-level, the PT-distribution was found to be sensitive both to the 
event geometry (the angles between jets) and the colour structure of the underlying 
hard event (specific hard parton configuration). 

In this paper we study the leading NP corrections to the above PT results. The 
NP contributions originate from the part of the phase space that corresponds to 
emission of secondary gluons from "large" distances, where the QCD coupling runs 
into trouble and the PT treatment becomes questionable. Our present understanding 
of QCD does not provide a "fundamental" solution to the problem. For a CIS 
observable, however, this part of the phase space does not lead to divergences (neither 
coUinear nor soft), but gives rise to a genuine NP contribution. The latter cannot be 
calculated from "first principles". The only firm information we have is that long- 
distance contributions bring in corrections to the PT results that are suppressed as 
an inverse power of the hard scale Q. The exponent p of the suppression factor 
can be extracted from the analysis of the divergence of the formal PT series 
(renormalon analysis, for a recent review see [4]). 

To make quantitative predictions and to make it possible to compare NP contri- 
butions to different observables, it has been suggested [3] to assume that the running 
coupling can still be defined at large distances. This simple assumption allows one 
to give meaning to PT contributions at any order and to relate the magnitudes of 
the power corrections with momentum averages of the QCD coupling in the infrared 
region. 

We follow the dispersive method [8] which extends the notion of the QCD cou- 
pling to small momentum scales and helps to extract and quantify the large-distance 
contribution to a given CIS observable. We show that the leading NP correction to 
Xout distributions and means belongs to the same class as the two-jet shape observ- 
ables listed above [5]. 

The main results, both in PT and NP sectors, were announced and their physical 
aspects anticipated in a short letter [6]. In this paper we present the derivation of 
the NP contributions to K^^^ -distributions and means. 

The results for the ifout-distributions that are common to the 1— T, M'^ /Q'^, C 
and B distributions are: 

• the leading NP corrections are of order 1/Q. Since the actual expansion param- 
eter is l/iTout rather than they are leading as long as -K'out is larger than 
Aqcd- For i^out of the order of Aqcd, higher powers become equally important 

[9]; 

• the leading NP corrections in the region Aqcd ^ -^out ^ Q can be embodied 
in a shift 5Kovx of the argument of the PT distributions; 

• the shift is a logarithmic function of Xout and is determined analj^ically by 
perturbative calculations, except for the overall normalization parameter A'^^; 
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• the phenomenological NP parameter A is given by the momentum integral 
of the QCD couphng in the infrared region. It is the same parameter that 
determines NP corrections to the two-jet shape observables (T, M^, C, B). 

Specific features of the X^t^-distributions are the following: 

• for the right X^^-distribution, the shift behaves logarithmically as In^^, 
which behaviour is similar to that of the single-jet broadening [7] ; 

• for the total iT^^-distribution, the logarithmic contribution In is supple- 
mented with a singular constant term proportional to 1/ which contribution 
has a weird colour structure and dominates at moderately small ii^out- This be- 
haviour is due to an interplay between FT and NP effects and is reminiscent 
of the total broadening case [7] ; 

• the relevant hard scales Q^^ are affected both by the small- and large-angle 
radiation of secondary partons. For the K^^^ distribution, these scales have 
a simple geometrical interpretation. As a result, the NP corrections, as well 
as the PT distributions, depend on geometry and colour structure of the hard 
event. 

A certain similarity between the i^out and the broadening distribution is due to 
the fact that both these observables accumulate contributions from arbitrarily large 
soft parton rapidities (while T, M^,C are dominated by finite rapidities — large- 
angle radiation only). As found in [7], it is this feature that is responsible for the 
logarithmic behaviour and the singular contributions to the shift. 

The need to present reproduceable derivations can sometimes be detrimental to 
readability of a paper. This is especially true for the present topic where the tricky 
event-plane kinematics and the NLO accuracy lead to technical complications. We 
have nevertheless tried our best to keep the paper readable. 

In section 2 we set up the kinematics, recall the PT results that will be needed 
for the NP analysis, and describe the NP corrections to the Sudakov exponents 
(radiators). Sections 3 and 4 contain the new results of this study, namely the 
NP corrections to the distributions and means of K^^^ and K^^^ observables. In 
section 5 we summarize the results and discuss their physical properties. The detailed 
calculations are presented in a series of Appendices. 

2. Distributions and resummation 

In the near-to-planar three-jet region (1.1) the event, at parton level, can be treated 
as a hard quark-antiquark-gluon system accompanied by an ensemble of secondary 
partons ki. We denote by Pa {a — 1, 2, 3) the three energy-ordered hard parton 
momenta: pi > P2 > Ps- 
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Defining the event plane as the {yz}-p\ane (we set the thrust axis and the thrust- 
major axis equal to the z- and y-axis, respectively), the total out-of-plane momentum 
variable is defined by 

3 n 

kl = Ei^-i + Ei^-i • (2-1) 

a=l 1=1 

We consider also the observable restricted to the right region (R) , that is the hemi- 
sphere with the lowest broadening: 

= I Pi. I + E I I • (2.2) 
ieR 

The hardest parton pi lies in the right hemisphere, near the thrust axis. 

In this section we first set up the kinematics of partons in the region (1.1) 
and define the event plane. We then recall the PT results that are needed for the 
NP analysis, in particular the factorization of soft emission and the result of SL 
resummation (details are reported in Appendix A). Wc finally report and discuss the 
NP corrections to the Sudakov exponents (radiators) that are evaluated in Appendix 
B. 

2.1 Pcirton kinematics and the observable 

At Born level one considers only the hard qq, g system and neglects the soft radiation 
so that Kont — 0. We denote the three hard massless Born momenta by Pa, take Pi 
along the thrust axis and P2 more energetic than P3, see Fig. 1. Given T and Tm we 
have^ 

Pi = £;(r,o,o,r), P2 = £^(^2,o,rM,-t2) , P3 = ^(2^3, 0, -Tm, -^3) , 



2-T T T 2-T Tlj Q (2-3) 

^^^^+2^' ^^^2+^^' P^Y'T^^ ^^2- 

The fractions 2:3 and ^3 are obtained from 2 = T + X2 + Xs and T = t2 + t^. 

There are essentially three underlying configurations for the hard partons ac- 
cording to which is the parton momentum corresponding to the gluon. We denote 
by 5 the configuration in which the momentum of the hard gluon is P^. In Fig. 1 we 
represent the three Born configurations. 

Beyond the Born level, the three hard parton momenta Pa no longer lie in the 
plane but are displaced from P^ by "soft recoils" q^. 

Pa^Pa + qa- (2.4) 
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To define the massless momenta Pa one requires that T and Tm are restricted to the region 

2(1-T)V2T- 1 - 

^ <Tm<VT^. 

In the foUowing we consider T, Tm restricted to this region. 
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Figure 1: The three Born configurations 5 for T = 0.75 and Tm = 0.48 ordered according 

to increasing probability. The thrust T and thrust major Tm are along the z- and y- 
axis respectively. The up, down, left and right hemispheres {U, D, L, R respectively) are 
indicated. 



As we have shown in [1], finite rescaling of the in-plane momenta Pa^, P2y,P3y, due to 
hard coUinear splittings, gets absorbed into the first hard correction to the emission 
probability of soft gluons, which is then resummed and embodied into the radiator. 
Bearing this in mind, all recoils and the secondary parton momenta ki can be 
treated as small (soft). 

Due to the fact that paz ~ Paz = 0{Q), p2y P2y = (Q) and p^y ~ Ps.^ = 0{Q), 
the five soft in-plane recoil components qaz,Q2y,Q3y can be neglected both in the 
definition of the event plane and in the soft matrix elements. Then, the definition 
of the event plane involves the remaining four recoil components and results in the 
following kinematical constraints: 



Qlt + ht — , Qit = {qix, Qly) , kit = {kixi hy) 

^ _ (2.5) 

q2x + qtx + ^'^ix = ^, (l3x + qix + ^f^ix = , qfx = Qlx ^{±Qly) , 



U D 



where we denote by R, L, U and D the right-, left-, up- and down- phase space 
regions, see Fig. 1. The total and the right-hemisphere Kout become 



i ieR 

2.2 Resummed distribution for K^ut ^ Q 

We study the integrated Xout-distribution defined as 



(2.6) 



dTdTM 



m \heR ) \heU J \ h=l 



(2.7) 
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where m denotes the number of final particles with momenta ph- The two delta- 
functions fix the event plane and the theta-function defines the observable. For the 
total- i^Tout case {Kj^^) the sum runs over all particles in the event. Analogously 
we define the right distribution by restricting the sum over particles in the 
theta-function of (2.7) to those belonging to the right hemisphere. 

In the near-to-planar "three-jet region" (1.1), the distributions in or in K^^^ 
are factorized in the form 



c[a. 



(2.8) 



= 1 

The factor daf^ /dTdTM is the three-jet differential Born cross section in the config- 
uration 6. The factor sj^^ accounts for the soft radiation emitted by the hard qq, g 
system. The first factor 

c{a,) = l + (a,) 

is the non-logarithmic coefficient function which takes into account corrections due 
to hard large-angle radiation not included in 'E'^^^{Ko^^l). The distribution sj*^^ is 
given by 

^[/^(i^out) = E ^ / ^^n/^(^out) n ^ Ks ■ (2-9) 

n I 

Here g is the distribution for the emission of n soft partons from the primary 
qq,g system in the configuration 5 and dHn (-K'out) is (a part of) the phase space, 
involving the three-jet kinematics (2.5) and the observable (2.6). 

The procedure used to resum the PT contributions of the Xout-distribution is 
described in detail in [1] and recalled in Appendix A. The essential point which 
allows resummation of 'E^^^{Kout) is the factorization of the soft distribution M^g 
and of the phase space dHn^^{Kout)- To achieve SL accuracy one needs to consider 
the soft emission factorization of M^ g at the two-loop level. 

Resummation employs the Mellin-Fourier representation of T/^^^{Kout) and gives 
rise to the exponential of the radiator. To describe this structure (see [1], Appendices 
A and B), we start from the phase space factor dH^ which is given by 



dH^(Kont) = d^qudq2xdq3x ■& l-f^out-E I l~E ' 

\ a=l i=l 

X S"^ ^it + ^ki^ s(^2x+qtx+^f^^^ (>(^3x+qix + ^kix^ , 



(2.10) 
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where the deha functions define the event plane (2.5). Correspondingly, the phase 
space factor dH^ involves only partons emitted in the right region: 

dH^{Ko^t) = dq^^ i& (x^^t -\qix\-Y.\''i^ ^ + j ■ (2-11) 

In the latter case, since the distribution is insensitive to the soft recoil components 
(l2x,(l3x and qiy, they can be neglected in the matrix element and integrated over. 

To factorize the multi-parton emissions in the phase space factors dHn^^ we use 
the MeUin and Fourier representations for the theta- and delta-functions in (2.10) 
and (2.11). This results in 

= / ^^^^ ' = / ^^^"^ n ^^"'^^^) ' (2-12) 

where the integration elements dji^^^ and the parton sources u^/^{ki) are described 
in Appendix A and B. Here we report some essential points concerning the integra- 
tion and the sources. 

Integration. For the Kj^^ case we need to introduce four Fourier variables 7, /?„ 
conjugated with qiy and qax, respectively. For the case only one Fourier variable 
P conjugated with qi^; is needed. The two integrations are 



8 ./_oo 27r ^J[J-oo 



(2.13) 

J ^ = 2 y — ~ — ^ ^1^) ■ 

Here we have rescaled the Fourier variables by v, so that they are dimensionless. The 
functions V'^^^ are given in (A. 7) and (A. 20). They are the Mellin-Fourier factors 
involving the recoil momenta. 

Sources. The soft-parton sources uF^^{ki) contain MeUin - Fourier factors involv- 
ing the secondary parton momenta. The sources enter as the factors [u{ki) — 1] 
integrated with the radiation probability dw{ki), where the subtraction unity corre- 
sponds to the virtual contribution. The source function tends to unity, u{ki) — > 1, 
when ki becomes collinear to one of the radiating hard partons Pa- This ensures 
cancellation of the singularities in the soft matrix element M^^. 

As has been already mentioned, for the sake of the PT-resummation the parton 
recoil effects can be neglected in the soft matrix element [7], so that we can take 
u{ki)^l for ki collinear to the skeleton-momenta Pa (see (A. 6) and (A. 21)). In the 
calculation of the NP corrections, however, it is essential to keep the recoil in the soft 
matrix element, so as described in detail in Appendix B. This makes it necessary 
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to accurately define the sources with account of the recoil as well, as to preserve the 
real- virtual cancellations in the collinear limit, see (B.19) and (B.34). 

The i^Tout-distributions share this feature with the broadenings [7]. In both cases 
the observables are uniform in secondary parton rapidities. In our case, accumulation 
of the NP contribution to Kout over large rapidities is bounded by the mismatch 
between the actual parton momentum and the corresponding "Born" momentum 
Pa (which lies in the event plane). Thus, the rapidity integrals turn out to be 
essentially the logarithm of the angle between the hard parton #a and the event 
plane. 

From (2.12) we obtain the Mellin representation of (-ft^out)- 

Sr(A„0 = /^.rM, (2.14) 

where the contour runs parallel to the imaginary axis with Reu > 0. The MeUin 
moment aj^^{i') is given by the sum, see (2.9), 



" ■ ' ' , (2.15) 



where the radiator TZ^/^ at SL level is obtained by using soft emission factorization 
of g at two loops. The radiator depends on the hard momenta Pa, the MeUin and 
Fourier variables and, in general, on the recoil momenta. 

The Mellin moment aj^^{u) is a regular function of u in the entire complex 
plane. It decreases for Reu — > oo so that 

Sr''(i^out<0) = 0. (2.16) 

The fact that al^^{v) cx e^^ont for R&v — > — oo, where K'^^ is the maximum kine- 
matical value of K^JJ^-, ensures that 

Sr^(i^out > KZ^ = <yV\^) = 1 • (2-17) 

The near-to-planar kinematics Xout <^ Q corresponds to the region of the Mellin 
variable Qv ^ 1. 

2.3 NP correction to the mean 

From the Mellin representation (2.14), we formally get 

- r KdK '^!^ = -a„a["'(.)|„„ , (2.18) 



9 



where, due to (2.17), the integration over K has been extended to infinity. This 
quantity, at the PT level, is determined by the region -ft'out ~ Q where the soft 
approximation is invalid. However, the NP correction originates, in any case, from 
the region of small (transverse) momenta of the order of Aqcd- Therefore the NP 
correction to the mean K^l^ can be computed within the present formalism based 
on the soft factorization. The NP correction to the mean is then obtained, according 
to (2.18), from the NP contribution to (T^'^iy) at i/ = 0. Even though Qv ^ O (1), 
the NP answer is still dominated by small secondary gluon momenta, that is, by 
large values of the Fourier variables /?„ and 7: Qv(3a ~ Qv^ ^ 1. 

2.4 -ftTjut case: NP radiator and recollection of PT results 

We discuss the structure of the PT and NP components of the radiator for the K]^^^ 
distribution. The detailed calculation of the NP corrections are given in Appendix B. 
We also recall the essential points of the PT results both for the radiator and the 
distribution (see [1] and Appendix A), which are needed to derive the NP corrections 
in the next section. In this subsection we suppress the index T and the index 5 
marking the jet configuration. 

2.4.1 Xj^t-radiator 

The radiator for the distribution has PT and NP components: 



The PT contribution depends on the Fourier variables but is independent of the 
recoil momenta. (As mentioned, recoil effects at PT level produce corrections that 
are beyond the SL accuracy, see [7]). On the contrary, the NP correction depends on 
the out-of-planc ,x-components of the recoil momenta but does not depend on 7, /3a. 

We first address the NP contribution to the radiator calculated in detail in 
Appendix B. We have 



We comment on this expression. 

• C(j is the colour charge of parton #a in the configuration 5 (here implicit). For 

instance, for the most probable configuration 5 = 3 in which the gluon is the 
least energetic of the three hard partons, we have C1 — C2 — CF and Cs — Nc- 

• Expression (2.20) is obtained by linearizing the sources u^{ki). The NP pa- 
rameter A^^ defined in (B.30), has dimension of mass and is responsible for 
1/Q power corrections to jet observables. 



7^(I/,A7,g)=7^PT(^^,A7) + ^^(?)• 



(2.19) 




(2.20) 
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• Linearization makes 51Z proportional to the Mellin variable f which, in turn, 
implies that the NP correction emerges essentially as a shift in /Tout- 

• The logarithmic behaviour in \qax\ comes from the soft gluon rapidity integra- 
tion which is bounded by the rapidity of the parton p^, with rapidities defined 
with respect to the corresponding Born momenta Pa- (Radiation at smaller an- 
gles corresponds to collinear splitting and does not contribute to the observable 
due to real-virtual cancellation.) 

• The hard NP scale Q^^ ~ Q associated with hard parton #a contains two 
factors: 

Qr = Qa•C^^ C^^ = 2e-^ (2.21) 

The rescaling factor C^^^ is defined in (B.27) and originates as a constant leftover 
after the real- virtual cancellation of collinear singularities in the NP radiation. 

The first factor Qa is determined by large-angle soft emissions and is sensitive 
to the geometry of the event as it depends on the angles between jets. In the 
configuration 5 one has 

= = 2(nP.) , Q] ^ ^^^'^p]ff'^ , (2.22) 

where b,c ^ 6 denote the quark- ant iquark indices. The scales for the fermions 
are equal to the invariant mass of the quark- ant iquark dipole. For the gluon, 
Qs is the invariant gluon transverse momentum with respect to the qq pair. 
(Notice that when the gluon becomes collinear to the quark or the antiquark, 
the scale Qs decreases and the non-Abehan contribution reduces.) 

• The NP radiator does not depend on the Fourier variables f3a, 7- As in the case 
of broadening, see [7], this is due to the fact that the contributions linear in 
kx or ky vanish, and the 7-dependence of the radiator emerges only at the 
1/Q2 level. 

• Since the NP radiator depends on qax, the way it affects the i^out-distribution 
depends on the PT distributions in parton recoils. 

The NP corrections to the integrated distribution I!,{Kout) that correspond to the 
NP radiator STZ are analyzed in the next section. 

Now we recall the structure of the PT radiator in (2.19) (for details see Appendix 

A): 

7^PT(l/, /3, 7) = Ri(P, /3i2, /3i3, 7) + R2{P2) + RsiPs) , 

r°° dv (2.23) 

Ri(^>,A2,A3,7)= / ^A Ri{Pi2) + Ri{pi^} 

Jo 7r(l+y^) 
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with = /?i + Pa) 



p^ = u^l + p2, a = 2,3, P^e^^u. 

The three terms of TZ^'^ in (2.23) describe the radiation from partons #1, #2 and 
#3, respectively. The two pieces of Ri correspond to secondary gluon emission in 
the Up {ky > 0) and Down hemispheres {ky < 0). The function Ra{p) is the two- loop 
single-jet radiator given by 

Raip) ^ Cavip, QD , r {p, QD = - [ ^ ^a,{2k^) In ^ , (2.25) 

with the running couphng in the physical ( "bremsstrahlung" ; CMW) scheme [10]. 
(The factor 2 in the argument of the running coupling originates after the integration 
over the in-plane component of the transverse momentum.) The perturbative hard 
scales are again composed of two factors: 

QV-Qa-Cr, Cr=b-'/'fora^S{qorq), ^ = le'^/^^^ for a = 5 ( ^ ). 

(2.26) 

The first are the same geometry- dependent factors (2.22) as in the NP scales, that 
account for soft large-angle emission. The second factors are hard SL corrections due 
to coUinear quark {—3Cp) and gluon (—/So) splittings. (The factor 1/2 in again, 
comes from the integration over the in-plane components.) 

The function r contains both DL (q;"L"+^) and SL (q;"L") contributions, with 
L — hi pQ and n > 1. The precise value of the hard scale and, in particular, the 
geometry of the event, is relevant only at SL level. 

The same is true for the dependence of the functions p in (2.24) on the Fourier 
variables f3a,'y, provided the latter are finite. This is the situation one encounters 
when calculating the PT distribution. Here the 7-integrals are well convergent 
so that /3a,7 = C(l) and we can simplify the 7-dependence by expanding the 
logarithmic functions Ra in (2.23) to SL accuracy as 

r{p, QD = r (P, QD + In (0 r'{u, QD + O (a.) , (2.27) 

where 

r'(i/, Q) = ud, r{v, Q) = ^^illM in(i.g) . (2.28) 

TT 

Since r' is a SL function, we do not need to care about the precise expression for the 
hard scale in (2.28) and about the difference between P, u/2 or ahke. 



12 



Physically, the dominant integration region /?a, 7 = C> {1) corresponds to the 
situation when the hard parton recoils and the secondary parton momenta are of the 
same order: 

Using (2.27) we derive then, with SL accuracy, 

-RF^iv, 7) ^ 5^ Ra{i^) + S{v, (3, 7) . (2.29) 

a 

where S is the SL function 

^ = r^i^irh^^ V((i-^7y)H/?L)((i+^7y)H/?i^3) 

""^^^^^ ' , (2.30) 

+ J]i?»lnv^, 

o=2,3 

and R'a{i^) — Car'{-u,Q). In the expression (2.29), the geometry of the event enters 
through the DL radiators Ra{v) while the SL function 5" contains the Fourier variable 
dependence and embodies hard parton recoil effects. Through the colour factors in 
(2.30), the recoil effects depend on the specific configuration 6. 

Let us stress that the approximate expression (2.29) can be used only if the 
/5a, 7-integrals converge, which is the case of the PT distribution. In the evaluation 
of the NP correction instead, one encounters poorly convergent (logarithmic) (3a, J- 
integrals and we cannot use the approximation (2.29) but have to employ the full 
expression (2.23). 

2.4.2 Recollection of the PT distribution in K^^^ 

Since the PT radiator is recoil-insensitive, we can freely integrate over qiy and Qax- 
The integrals in the Fourier variables are convergent, so that we can use the approx- 
imation (2.29) and obtain, to SL accuracy, the final result 

a^'^iiy)^! dA^^e-^''''^'^'^'^) = 6"^ ^»(^) JS^(i/) , Tt{v)^ j dfj,'^ e-^^"^^^"^ , (2.31) 

where dfi^ is the integration introduced in (2.13). Here (Ft{^) depends on v via the 
SL function r'(i/, Q) and has been computed in [1]. To first order in r' one has 

J^t{i^) = 1 - In 2 r' ( 2Ci + C2 + C3 ) + O (r") R!^ = C„ r'{v, Q) . (2.32) 

We now recall how to compute the distribution, i.e. how to perform the i/-integral 
in the MeUin representation (2.14). To this end we use the identity 



r(i + z) 



z=0 

(2.33) 
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The action of this differential operator on the function r~^(l + z) reduces, within 
the adopted accuracy, simply to replacing z with Rip. The PT distribution to SL 
accuracy takes the form 

E-^(ir,„0 ^ e- S» «"(^o-^t) . , R'^^Cr r' (i^^ , Q) , (2.34) 

with Ct — 2Cf + Nc the total colour charge. Dependence on the geometry enters 
through the hard scales in the DL radiators Ra{Kout): while the hard parton recoils 
build up the SL correction factor J^TiK~^^). 

The role of recoil effects and the event-plane kinematics can be seen already in 
the first-order expansion in as. Using the expansion of J-t given in (2.32), for the 
differential distribution we obtain 

—nr^ = — i^il^:^ + ^2ln-^ + C3ln-^^+ 0{a,) . (2.35) 

Gluon emission off the parton has its characteristic hard scale Qa^. The additional 
factors 4 and 2, that enter the hard scales, originate from the SL function J^t in 
(2.32) and have a simple kinematical origin. The standard definition of the event 
plane implies that when the secondary gluon is emitted in the right hemisphere (see 
Fig. 2b), all three hard partons experience equal recoils, 

k,>0, ky>0: pi^ = -k^ = p2x = -psx \ , , , K^^t 

ky <0: pi^ = -k^ = -p2x =P3x j 4 

On the other hand, for the secondary gluon in the left hemisphere (see Fig. 2c and 
2d), only one parton recoils against it: 

k^<0, ky>0: p2x = -kx] Pix = Psx = \ K^nt ,^ 

( =^ \kx\ — ■ (2.37) 
ky<^: pso; = -kx\ Plx =P2x = j 2 

These are the factors 4 or 2 entering the scales in (2.35). 

2.5 i^(^ut case: NP radiator and recollection of PT results 



Here we describe how the results of the previous subsection change in the case of the 
distribution. (In this subsection the suppressed indices R and 5 are implied.) 

2.5.1 X^„t-radiator 

As in the previous case, the radiator for the K^^^ distribution has PT and NP com- 
ponents: 

n{u, (3, gi,) = n^^iu, (3) + 5n{qix) . (2.38) 
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(e) (d) 

Figure 2: (a) Hard partons in a generic Born configuration, (b) Soft gluon k (the short 
curly stick) is emitted in the right hemisphere. All three hard partons experience equal 
out-of-plane recoils, see (2.36). White (shadowed) partons have positive (negative) x- 
components of the momentum, (c) Soft gluon k is emitted in the up-left region. According 
to (2.37), only parton ^2 recoils (shadowed; it has a negative j;-components of the mo- 
mentum), (d) The case of k in the down- left region. 



The PT radiator does not depend on gi^. while the NP correction does not depend 
on p. 

The NP correction is described in details in Appendix B. Since the observable 
accumulates only particles emitted in the right hemisphere which contains only one 
hard parton #1, the NP radiator can be written as 

-oNP 

6n{qi,) = z/A^^Cilnf^. (2.39) 

l^lxl 

Subleading contributions due to radiation into the R hemisphere off the L-hemisphere 
partons are embodied in the hard NP scale Q^^ ~ Q. This new scale is not related 
to the NP hard scale (2.21) for the K^^^ case.^ Contrary to the K'^^^-case, the hard 
scale in (2.32) does not have a simple geometrical interpretation. This is due to the 
fact that particle selection in the i?-hemisphere is "unnatural" for the two-parton 
antenna dipoles that determine the structure of the accompanying soft radiation in 

^Denoting, for the sake of simplicity, the hard scales for R- and T-distributions by the same 
symbol does not cause confusion. 
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the three-jet event. The procedure for computing the i?-scale in (2.39) is described 
in Appendix B.5. 

To evaluate the effect of the NP radiator upon the distribution we need to know 
the PT radiator. It has been computed in [1] (see also Appendix A) and reads 

W^iuP) = R^(p) = Ci r(p, QD , P = Pa/TT^ , (2.40) 

with r the DL function defined in (2.25). The hard PT scale ~ Q associated with 
the soft parton emission in the right region, does not coincide with the hard scale 
(2.26) for the PT i^T^^-radiator and, once again, does not have a simple geometrical 
meaning. 

2.5.2 Recollection of the PT distribution in K^^^. 

The situation here is simpler than before since we have a single Fourier variable f3. 
Neglecting the NP radiator, the gia;-integral can be freely performed and one obtains 

By using the identity 

/W^^,Jfg£L. (2.42) 

we obtain the SL expression of the distribution 

E-^(ir„„0 c e-MK^) . Y§^y R'l ^ (i^- , Q) , (2.43) 
with corrections beyond SL accuracy. In general we introduce the jet #a function 

/■oo JO -Ra{^>^/^) 

W = e^"(^) / ^° ^ f{K{^)) . R'M ^ Cy{u, Q) , (2.44) 

As before, the O (ctg) expansion of the differential distribution shows that the SL 
corrections have a simple kinematical explanation. Prom f{R'a) — 1 — ln2 • we 
obtain 

^^P^A,,J-^+0(ar). (2.45) 

The additional factor 2 in the scale has the following interpretation. When a single 
gluon k is emitted in the i?-hemisphere (see Fig. 2b) all four partons, as we already 
know, experience equal recoils, \kx\- Counting only the i?-hemisphere particles we 
have Kout = 2|A;^|. 
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3. NP correction to the K^^^ distribution and mean 

(In this section the suppressed indices T and 5 are imphcit, unless exphcitly needed.) 
Taking into account the NP radiator (2.20) we obtain, in the hnear approximation 
in l/Q, (see Appendix B) 

I Q.a,x I 



a 

= j d/.^e-^"^('^'/^'-)|l-^A^^j:a[ln(pgf)+X„(/3,7)] 



(3.1) 



The unity in the curly brackets gives the PT contribution, while the second term 
describes the leading NP correction. The functions Xa ai'c given in Appendix B. 
They depend only on the Fourier variables and are constructed with use of the 
function x(/3): 

x(/3)=lnv^lT;3^+/3tan-i/3=|/3|^+x'(/3), H x\P)-0, (3.2) 

with x'iP) = large f3. For a = 2, 3 we have Xa = xiPa)- The function 

Xi(/3,7) contains the combination of xiPu) and xiPis) and is defined in (B.40). 

Equation (3.2) shows that the functions linearly grow for large \ Pa\, and 
this makes the structure of the NP corrections quite different from that of the PT 
contributions. Indeed, for the PT contribution the Fourier variable integrations in 
d/ji^ are rapidly convergent, so that all soft momenta are of order Kout- For the NP 
corrections, due to the presence of the Xa functions, the Fourier variable integrations 
in dfi^ become logarithmic. This implies that the phase space regions with the 
secondary gluon momenta much smaller than i^out become important. 

From (2.31) and (3.1) we obtain 

a{u) = a^'^iu) - 5] C„ fa{u) , (3.3) 

a 

where 

Uu)^a''^{u) ■ In(pg^P) + Uu) , Uu) = j df,^ UP, t) e'^''^^^''^'^) . (3.4) 

The functions /„ are analyzed in Appendix C. We show that they are given by an 
integral similar to the one for the PT contribution to the Mellin moment, except for 
the presence of the x functions in the integrand: 

3 



J —oo '^[^ ~r Pa) 
^-oo ~r Pa) 



(3.5) 
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where (323 = (32 + (33- The factor S is the contribution of the parton #1 (see Ap- 
pendix A) and is given by 

S{u, /3, P23) = e-^^(^'^'^-°) -I rtLi^ g-Ri(pA/3.,7) ^ = p+p^^ ^ (3.6) 

Jo 7 

with Ri defined in (2.23). 

The presence of the function xiPa), hnearly growing with \(3a\, in the integrand 
of (3.5), makes the evaluation of fa{^) quite involved. The large values of Pa, J, 
corresponding to the recoil momenta much smaller than Kout; give the dominant 
contribution to /„. In these circumstances we cannot use the SL expansion of the 
PT radiator (2.29) which simphfied the calculation for finite values of the Fourier 
variables. 

The strategy of the calculation is the following. We split xiPa) into two pieces 
(3.2). The contribution of x'iPa) is evaluated by the same techniques that we used 
for the PT case. Indeed, since x'iPa) is growing only logarithmically, the Fourier 
variable integrations remain convergent and we can use the SL expansion (2.29). 
The linear piece ||/3a| which gives the leading contribution to /a(f^) will be analyzed 
separately. 

Result for /i. In this case the /3i-intcgration is logarithmic while the remaining 
1,(^2,^3 integrals remain convergent. To analyze the leading piece we need to study 
the behaviour of S for large Pi and finite ■j,P2,P3. In this limit the imaginary part 
of Ri vanishes, and S reduces to the first term only, see Appendix C.l, 

S{u,Pi,Pt) ^ e-«^(^^) , pu, Pi3 -Pi = i^^l + Pl (3.7) 

This behaviour can be simply understood by considering the event-plane kinematics 
(2.36), (2.37). For /^i — > cxd we are forcing parton #1 to stay in the plane by 
suppressing PT radiation off it. At the same time, radiation off the other two partons 
is not constrained, since parton 7^1 does not recoil when PT emission occurs in the 
L- hemisphere. As a result, S turns out to be simply the Sudakov suppression factor 
for not having PT gluon radiation off the parton #1. 

The leading piece of /i is obtained by substituting xiPi) ~^ f I'^il ("^■^) 
using the approximation (3.7). All Pa integrals are now factorized and we obtain 



a=2,3 



ne-^"(^^^2(^)^3(i^)£^i(f?), 



a=l 



with ^a(i^) the SL functions defined in (2.44). The function Ei{u) given by 



E,{u) = ^ e-«^(^)+^^(^) (3.9) 
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has been studied in [7]. Its properties are recalled in Appendix D. 

Corrections to the leading behaviour of /i are discussed in Appendix C. The SL 
result reads 

/i(i.) = e-s«^»(^){:^2H^3M£^i(p) + CiH} , (3.10) 

with Ci(z/) a SL function which vanishes for z/ = and is proportional to r'(u, Q). 

The function Ei{K~^^) is related to the restricted average of InQf^/lo'i^l, under 
the condition \qix\ < K^ut- To compute this average, we use the Sudakov distribution 
in the recoil momentum q = \qix\ of parton #1, 

V,{q) = e-^^^^-'\ (3.11) 



and derive 



\ 191^1/ ki,|<i^out ^ 
r,NP 

In ^ + E,{K-:,) . 



(3.12) 



From the behaviour of the function Ei, see (D.4) and (D.5), we obtain the leading 
terms of Li 



L,{K,,) = ]n$^{l + 0{ (a^ln^-^) | 1, a^ln' Q / K,,, :^ 1 





(3.13) 



^i(i^out) = <^ln [l + O (a,\n' ^ ) ) , ajn' Q/K,^, « 1 

The unrestricted average is given by 



0_ 



^^th) = Ei[^]=-j===-\n^-^ + 0{^,) . (3.14) 



Up to corrections of order ^/a^ it is independent of the event geometry. 

The restricted average Li{Kout) is a decreasing function of -fTout- It starts from 
InQ^^/Kont for very small K^ut (^s In^ Q/Xout > !)■ With K^ut increasing, it 
freezes at the value of the unrestricted average Ei{l/Q^^) at moderately small Xout 
(agln^QZ-fCout <S 1), up to corrections of O {^/oi^). 

Result for /„ with a = 2,3. Consider for example the case of f2- Here the 
P2 integration is logarithmic while the /?i,/?3 integrals converge. Evaluation of the 
leading behaviour of S for /^a <S P2 is subtle. In this /3-region the second term of S 
in (3.6) is not negligible any more since the 7-integration turns out to be logarithmic. 
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As shown in Appendix C, the leading exponential contributions exp[ — ^Ri) in the 
two terms of S cancel, and one is left with 

S{u,P,P,s)^e-<'^^. (3.15) 

This behaviour has the following physical explanation. For ^ oo we force p2 to 
stay in the plane. From (2.36) we see that in this case also pi has to be devoid of 
PT radiation, since gluon emission off parton #1 would cause parton #2 recoiling. 
Therefore, the distribution for large /?2 is given by the product of the two Sudakov 
form factors: one from parton #2, giving exp(— i?2(i^-\/l + /3|)) and the other one 
from the parton #1 giving (3.15). Radiation off the parton #3 is not constrained. 

The leading piece of /2 is obtained by substituting in (3.5) xiP2) f |/?2| and 
using the approximation (3.15). All f3a integrals in (3.5) are now factorized. The 
/9i-integral gives 1; the remaining two integrals produce 



3 
a=l 

where we used p — V-^/l + /?|. The function E2 is defined as 

Ea{v)= /" ^ e-«i('')-^«('')e^i(^)+-^"(^) (3.17) 
Jd P 

with a — 2 and is analyzed in Appendix D.2. Corrections to the leading behaviour 
of /2 are discussed in Appendix C. The answer to SL accuracy reads 

h{y)=e-^''-^'^{:F,{u)E2{v)+C2{y)] , (3.18) 

with C2{v) a SL function which vanishes for z/ = and is proportional to r'{u,Q). 
The expression for is obtained from (3.18) by interchanging the labels #2 and ^3. 

The functions Ea{K^^^) with a = 2, 3 can be related, as above, with the average 
of ln{Q^^ /\qax\) under the condition \qax\ < Kont- 



g --Ri(p)g-R2(p) 



(3.16) 



\ / |?a.|</fout Jo Q Q 



^ln^ + E^{K;^,), 



(3.19) 



where the PT distribution over the parton recoil \qax\ is given by (the logarithmic 
derivative of) the product of the Sudakov form factors of partons #1 and #a: 

Vaiq) = -jf- e-^i(«")-^»(«") . (3.20) 
a In 5 
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Prom the behaviour of the functions Ea (see (D.7) and (D.8)) we have 

^a(i^out) = ^In (^1 + O (^asln' ) ) ' « ^ ' 

where the unrestricted averages now read 



(3.21) 



in;^ =E, 



a 1/ VONP 



(3.22) 



where we dropped corrections of absolute order (relative corrections 0{as)). 
The diagonal ratios Qa^/Qa^ include only hard corrections to the scales. The non- 
diagonal ratios Qi^/Qa^j the contrary, depend also on the event geometry. 

3.1 Correction to the mean K^^^ 

The NP correction to the mean is obtained from (2.18) and (3.3) by taking formally 
the hmit i/ ^ of the functions fa{^)-, which actually means vQ = 0{1). We get 



where we used the fact that, for v = 0{1/Q), we have Ra{i^) — ^ 0, J-'a{i^) — ^ 1, 
Ca(z/) and 1. The quantity in the curly brackets is Laii' ^). We 

obtain, up to corrections of order ^/a^, 

/„(0)^/ln^\, (3.24) 

\ \Qax \ I 

with the average given in (3.14) and (3.22). The final result for the mean is then 

/ nNP 

\NP xNPY^A^ /i„'^a 



\ \(lax\l 

By using the expansions (3.14) and (3.22) we get 

(^o'^t)'''' = A^^ ^ f^/cT + , + ^^=1 + O (1) . (3.26) 

The legitimate constant term here can be restored from (3.14) and (3.22); it contains 
the ratios of hard PT and NP scales and depends on the geometry of the event. 

The peculiar colour structure of (3.26) shows that, due to an interplay between 
the PT and NP effects, the three jets do not contribute independently to the mean. 
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3.2 The distribution 

From the MeUin moment (3.3) one obtains the NP correction to the integrated dis- 
tribution 

S(i^out) = SPT(i^out) + 5S(irout) (3.27) 

where 

_^e'^^out/^(^,), (3.28) 

with fa{i^) given in (3.4). To evaluate the Mellin integral we use the operator identity 
= I d/.^e-^"^(-^^'5'-) jln g + ^(1 + .) - 1 + UP: 7)} 



r(i + ^) 



(3.29) 



z=0 



where in the first line we have replaced Inz/ with —dz- 

Consider the two pieces of ipa given by the first two and the last two terms in 
the curly brackets on the second line of (3.29). The first piece can be expressed in 
terms of the PT distribution as follows: 

V'i^)(Xo,0 ^ A« • 9^„,,EP^(X„,0 , (3.30) 

where 

^NP -qNP 

Ai^) = In 1^ + ^'(l + B^) = In 1^ (1 + O (a,)) . (3.31) 
Here we used the operator identity (2.33) in the form 

S^^(Xout) = / ci/.^e-^"^^^(-^^/^'-) , (3.32) 



z=0 



and substituted R'j, for z with corrections of order O (as). The contribution (3.30) 
to the NP correction 5E(ii'out) clearly has the form of a shift of the PT distribution. 

In a similar way we calculate the second contribution to 1pa{Ko^^l) given by the 
last two terms in the curly brackets of (3.29). Neglecting corrections of order ctg we 
have^ 

^(^)(Xo.O c -1- a^„„,E^^(Xo.o + ^ r° (/+ f^'^ • ^'-''^ 



^The substitution of z with i?^ is in general correct only for regular functions of z. Actually 
here we can apply the substitution for the function 1/z since the singularity cancels. With more 
cumbersome algebra one can avoid this intermediate step altogether. 



22 



Let us show that the second term with the derivative of /„ generates a l/i?^ contri- 
bution which cancels the singularity at R!j, = of the first term in (3.33). To verify 
this cancellation we consider the singular pieces of fa containing the Ea functions 
(see (3.8) and (3.16)). 

We start from the leading piece of /i in (3.8) and take the derivative to obtain 

e-Ei?.(i?-ut) 

Ok^J, ^ :F, {K-^,):F, {K-^,) {1 + {I^, + I^,)E, {R-^,) } . (3.34) 

-"■out 

Using (2.34) and (2.44), we arrive at 



r(i + R'^) j^j 



with all functions in the r.h.s. evaluated at K~^i-. The second term is singular for 
R'rp and cancels with the —l/Rip term in (3.33). (Notice that and tend 
to unity in this limit.) 

Similar results are obtained from the leading pieces of /2 and /a. Taking the 
derivative of (3.16) we obtain 

^i^out A (^o-ul) ^, J'S (C. , 1 



r(i + i?^) j^T [o 

and again the —1/Rip singularity in (3.33) cancels. 

The derivative of the subleading pieces Ca in the expressions (C.8) and (C.24) 
for fa is simpler to evaluate. We have 



dK^,. S^T-^ - ^5i.„,,S-^(i^out) , (3.37) 



'C p-T,bRb\ c 

w \^r(i + i?^) ) ~ Tt ^""^ 

where all functions are evaluated at K~^^ and we used again (2.34). 

In conclusion, the second piece i)if^{Kont) of the NP corrections of the K'^^^- 
distribution can also be written as a shift of the PT distribution. Up to corrections 
of relative order ctg, we have 

V'f (i^out) ^ Ai^) • a^„„,EP^(i^o,t) , (3.38) 



where 



A (2) C2 + Ca ^ 2T 3 T 2T 3 — Tt , C\ 

i^V ^ —7^ + y p. + 

c.^E, + ^^ + ^, (3.39) 

(2) C2J-'2 p T2—J^T C3 

3 — ^ -r- -^3 H ^ p/ ^ ^ ' 
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with all functions evaluated at -ft'out- The first terms are the leading pieces. The 
second terms are SL functions regular for i?^ (the numerators vanish when 
r' — >■ 0). The last terms are SL corrections which vanish with r' — > 0. 

Summing up the two corrections (3.30) and (3.38), Aa = aI^-* + A'^\ we finally 
obtain the expression for the NP shift 

E(Xo„0 = SPT(x„,,) + 5K,^t ■ 9^„„t^^^(irout) - ^""^{K^t+SK^t) , (3.40) 
with 

5Xout = -A^P J2 ^a(^out) . (3.41) 

a 

By expanding aI^"* and aI^"* and using the functions La{Ko^i) given in (3.12), (3.19) 
and describing the restricted averages, we arrive at 

Ai(i^out) ^ — Li{Kont) + 7^ In ^7 — , 

A,(ir„„0 ^ ^ L,(Ko.O + ^ In + ^ In ^ , (3.42) 

A3(is:out) - HKout) + 7^ In + 7^ In -r^ . 

-"-out <-^r -"-out 

Here we introduced the function La{Kout) = -^a(-f^out) — 7b and we neglected correc- 
tions of relative order r' = O (dg InQ/i^out)- 

To understand the structure of the result we consider two regimes. 

Qfs In^ QZ-ft'out 3> 1- This regime corresponds to well-developed PT cascades. From 
(3.13) and (3.21) we have La{Kout) ~ ln((5/-ft'out); so that the colour factors in (3.42) 
cancel and all Aa{K out) ~ ^T^{Q/Kout)- For the shift we then simply have 

SKout ^ -A^P E ^» + ^ (1)) ■ (3-43) 

In this regime multiple PT radiation causes all three hard partons to experience 
comparable recoil \qax\ ~ -f^out- The NP radiation off each of the partons then 
averages to InQ^^ /\qax\ — InQ^^/Xout thus producing the above result. 

as\n^ Q / Kout ^1. In this regime there are only few secondary PT partons. Let 
us restrict ourselves to the first-order in ag, that is to emission of a single PT gluon: 

E - c SKout ^ - 5i^out P- ^ 5i^out ^ . (3.44) 

"-"■out -"-out -"-out 
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Observing that i?y oc Ct, we can represent the characteristic product 6Kout R't 
the following equivalent form: 

3 

5Xout = E <(^o"ut) 8K^1 , (3.45) 
0=1 

where (up to ttg corrections) 

iif W = -A- |C: k ^ + C, fi^ + C3 In 

L -"-out -fiout -fi-out 

= -A^P |ci Li(Xout) + C2 In ^ + C3 hiK,,,)] , (3.46) 

= -A""^ |ci Zi(Xout) + C2 Z2(Xout) + C3 In ^ 

Being accompanied by the factor R'^{K~^^^_) oc Ca in (3.45), the contribution SK^^^ to 
the NP shift corresponds to the PT gluon emission off the parton ^a. 

Thus, the shift 6K^]^^^ given by the first line in (3.46) arises when the PT gluon 
is emitted from pi. In this case, as we know from (2.36), all partons undergo 
equal recoils, so that \qax\ ~ -^out/4. Each hard parton then contributes with 

The shift 5K^^^ arises when the PT gluon is emitted from p2- In this case 
only parton p2 recoils with \q2x\ — -f^out/2 (see (2.37)). In this situation, the 
NP gluon emitted from p2 contributes to the shift, as before, C2 In Q^^/|g2a;| — 
C2ln(2(5^^)/i^out (cf- the second term on the second line in (3.46)). When the NP 
radiation occurs off 7^1 and ^^3 the situation is entirely different. Since the partons 
^1 and 7^3 do not recoil against the PT gluon emission in the region #2 (up-left 
quadrant), the momenta pi and ^3 stay in the event plane. In these circumstances 
the NP radiation from #1 and #3 would diverge if not for the high-order Sudakov 
form factor effects that push the partons off the plane by certain 0. The cor- 
responding averages of InQ^^ /\qax\, a = 1,3 give the singular L-factors (3.12) and 
(3.19). In a similar way one can interpret various terms in SK^^. 

We notice in conclusion that, using (3.40), the X^^-distribution with account of 
the leading NP correction can be written in the form 



(3.47) 



where, up to correction of order r', 5i^out given in (3.46). 

4. NP correction for the K^^^ distribution and mean 

Here and in the rest of this section the indices R and 5 are implied. The integration 
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over the recoil qix gives (see (B.39) and (B.40)) 

r \ f ^ R -Ri(i>y/u^) f\Qix\\ 

+ {l-.A--C.[.n(.gn^xW]}. 

with corrections of the second order in A^^. From (2.15) and (4.1) we obtain the 
Melhn moment 

a{v)^a''^{v)-v\''''C,fM. /i(i^) = a^V) • MpQ^ + h{y) . (4-2) 
where j\ is given by 

= £ ^frf) '""'("^ = {^.W ■ (4.3) 

This is obtained by sphtting the x(/3) function (see (3.2)). The function Ei{v) is 
given in (3.9) (with the PT scale QX^ of the X^^. case), and Ci{v) is the SL function 

Here wc used the expansion Ri{V^/l+^) — Ri{iy) — In yT+^Ci r'{u, Q) since the 
/^-integration is fastly convergent. We now analyze the NP corrections to the mean 
and the distribution. 

Correction to the mean K^^^. The calculation is similar to that for the previous 
case. Prom (2.18), (4.2) and (4.3) we obtain 

{kL)T - ^"^^ MO) c C^ (in ^\ , (4.5) 



with corrections of order ^^a^. The average value here is computed over the parton 
#1 distribution given in (3.11) (apart from the hard scale being the one for the K^^^ 
case) . 

Correction to the K^^^ distribution. As before, from the Mellin moment (4.2) 
we obtain the NP correction to the distribution 

5E(Xout) = -A^P C„ Vi(i^out) , V'i(i^out) = dK,^^ J ^e'^^out /^(^) . (4.5) 
To evaluate ipi we use the operator identity 



(4.7) 

2=0 
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to get 



'dp / 1 



'2/ n^P 



out 



(4i 



z=0 



The integration over x(/3) generates a singular \/z term: 



J — ( 



dPxiP) 



TT 



l + Z 



dp 



TT 



1+2^ 



(4.9) 



This means that inside the integrand of (4.8) we can replace 



1 

h 

z 



m ^ (i + f ) - (^) 



(4.10) 



We conclude that the singular 1/z term in (4.8) is cancelled by the contribution from 
the integral of the x{P) function in (4.9). 

Performing the /9-integral in (4.8) we obtain the PT distribution {z is substituted 
with R[). As for the K^^^ case, the NP correction can be expressed as a shift 



(4.11) 



where 



Ai(Xout) = ln|i-+V^(l + i?;) + i^ (l + f ( 

(1 + C(a,)). 



= ln 



2Qr 
K. 



(4.12) 



out 



We see that in this case we do not have any contribution from the singular Ei- 
function. This is due to the fact that in (4.10) the 1/z singularity was exactly 
cancelled by the integration over x-function. In the K^^^ case the corresponding 
cancellation was only partial. 

The physical reason for this cancellation can be understood as follows. For 
the distribution we are considering emissions only in the right hemisphere. In 
this case the parton #1 always recoils with \q-ix\ ~ -f^out; be it the regime of well- 
developed PT cascades, as In^ Q/K^Mt ^ 1, or the regime with few secondary partons, 
ccs In^ (5/-f^out -C 1. (In this latter case even the finite factor 2 in the scale can be 
explained from the event plane kinematics (2.36) and Fig. 2.) 



5. Summary of relevant formulae and discussion 

In this paper we have computed the leading NP correction to the out-of-plane mo- 
mentum distribution defined in (2.7) in the kinematical region (1.1). We have studied 
the distributions in the variables (see (2.1)) and K^^^ (see (2.2)). The NP results 
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obtained here, together with the SL-resummed PT calculations carried out in [1], ex- 
tend the accuracy of QCD description of CIS jet observables in e'^e~ annihilation 
processes beyond two-jet physics. 

In spite of the relative complexity of three-jet ensembles, both PT and NP ex- 
pressions have a clear physical interpretation based on simple QCD considerations 
and kinematical relations. One should expect that the present methods can be ex- 
tended to other processes and CIS distributions in multi-jet environment. 

We have analysed the NP corrections to the logarithmic "soft factor" E(i^^out) in 
the general QCD expression (2.8) based on the QCD factorization of soft radiation. 
The corresponding non-logarithmic coefficient function is computed in [11]. 

The NP corrections have been obtained by using the dispersive method of Ref. [8] 
for defining the running coupling and triggering the NP effects in CIS observables. 
The characteristic NP parameter A'^^ defined in (B.30) can be related, after merg- 
ing of PT and NP contributions, with the integral of the running coupling over the 
infrared region. It is the same parameter that enters the NP corrections to distribu- 
tions and means of two-jet observables such as thrust T, invariant jet masses M^, the 
C-parameter and jet broadenings Bt, By/. In (B.32) A'^^ is related to the well-known 
NP shift in the 1— T distribution, see [5]. 

The previous sections are full of formulae, so it may be helpful to single out the 
important ones and to list their main physical properties. Before turning to the NP 
corrections, we recall, for the sake of completeness, the key features of SL-resummed 
PT distribution i:'^^{K^^^. 

5.1 PT results 

We discuss first the X^^. case and then K^^^. 

K^^^^ case. The PT result is given in (2.34) with the following structures: 

• the three radiators Ra{K~^i) in (2.25) collect all DL « In'^^^ Q/Kut) contribu- 
tions and essential SL corrections (a^ln" QZ-f^out)- Ra is the standard radiator 

that one encounters in the two-jet physics [2], slightly modified here by the 
factor 2 in the running coupling scale, which is due to the fact that the relative 
transverse momentum in the argument of the coupling is integrated over one 
(in-plane) component; 

• the hard scales Q^^ are given in (2.22) and (2.26). Their precise expressions 
are essential only in the arguments of the DL terms (the SL correction terms do 
not distinguish between Q and Q^^). It is important that these scales depend 
on the geometry of the event (the angles between jets); 

• the additional "recoil" factor J^riKout) is a SL function (a function of r' oc 
aslD.{Q / Kout)) given in [1] as a multiple integral over Fourier variables. This 
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function depends on the configuration 6 of the underlying hard process (Born 
parton system), as can be seen from its first-order expansion (2.32). The latter 
expression has a simple physical interpretation and can be easily reconstructed 
from event-plane kinematics with a single secondary gluon (see the discussion 
in subsection 2.4.2). 

K^^^ case. The PT result is given in (2.42) and (2.43). Here only one radiator 
contributes, Ri{K^-^^) of the most energetic parton in the right hemisphere, so that 
the distribution in K^^^^ is sensitive to the Born configuration already at DL level. 
Radiation from the other two hard parton contributes only at SL level and can 
therefore be included in the hard momentum scale Qi^, which has been computed 
in [1] and does not have a simple geometrical interpretation (see the argument given 
in subsection 2.5.1). As in the K^^^ case, the first-order expression for the distribution 
can be reconstructed by simply examining event-plane kinematics. 

5.2 NP result 

We have studied the NP contribution leading in A^^. This approximation is vahd in 
the region Aqcd <^ -f^out -C Q. For values of Xout of the order of Aqcd the simple 
shift approximation fails, see [9]. 

We start again from the case of the Kj^^ distribution. 

K^y^^ case. The distribution T,{Kout) with account of the leading NP correction, 
is given in (3.40). The NP correction enters as a shift of the argument of the PT 
distribution, SK^ut given in (3.41). The NP gluon emitted by the hard parton #a 
contributes to the shift with two terms A^\Ko^i) and Ali\Ko^t) given in (3.31) and 
(3.39) respectively. Corrections to these expressions are of the relative order ctg and 
are beyond our accuracy. 

The first term Aii\Kout) is additive (independent emission) since its contribu- 
tion to SKoat is proportional to the colour charge Ca of the emitter. The functions 
A''a\Kout) behave logarithmically in Kout, with the hard NP scales Q^^ given in 
(2.21) and (2.22). These scales differ from the corresponding PT scales (2.26) only 
by constant factors {(^^ versus (^'^) , which are due to small-angle radiation effects. 

The second term Aa (Kont) is expressed in terms of the PT functions Tt{K^^-^^), 
Ta{K~^^) and Ca{K^^^), and the NP function Ea{K^^^). In particular, Tt{K~^^ 
and Ta{K~^^) are the SL functions encountered in the PT analysis. The first is 
given in [1] and the second in (2.44). The SL functions Ca{K^^^) are defined in 
Appendix C. Finally, the NP functions Ea{K~^^)^ defined in (3.9) and (3.17), are 
typical for observables that are uniform in the rapidity of the emitted particles. 
They have first emerged in the study on the NP correction to the total broadening 
distributions in [7]. 

One has the following two regimes for the shift: 
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• for very small Kout (cis In^ Q/-^out ^ 1), when the QCD cascade is well- 
developed and there are many secondary partons around, the shift 6Kout is 
dominated by aI^-* , is logarithmic in Kont and proportional to the total colour 
charge Ct = 2Cf + Nc of the three-parton system: 

5Ko,t^-A^^CTln^. (5.1) 

In this region the NP radiation effects are additive; 

• in the opposite region with moderately small i^out (o^sln^ QZ-ft^out ^ 1); when 
there are few secondary partons, the shift is dominated by the contributions 

(2) 

Aa that are singular in ccg, and the dependence on the colour charges becomes 
rather odd: 

is- ~ JrA"'- f C,(C, + C,) , C,C, , C^C, \ 

This peculiar colour structure results from an interplay between the NP and PT 
radiation effects. Its physical origin is discussed in detail in subsection 3.2, below 
(3.42). 

K^^^^ case. The -ft^j^t distribution, with account of the leading NP correction, is 
given in (4.11) with the NP shift given in (4.12). The procedure for computing 
the hard NP scale Qf^ is presented in Appendix B. This scale does not have a 
simple geometrical meaning. In agreement with the usual argument based on the 
event plane kinematics (see (2.36) and the discussion after (4.12)), the shift in the 
distribution is purely logarithmic in K^vx- Unlike the i^^J^^ case, it does not 
contain singular 1 / ^/a^ contributions. This feature is similar to that of the single-jet 
broadening distribution [7]. 

NP corrections to the -fCout means. The results are given in (3.25) and (4.5) 
for the total and the right-hemisphere observables, respectively. They are expressed 
in terms of the unrestricted averages (3.14) and (3.22). Again, the peculiar colour- 
charge dependence is due to the PT-NP interplay. 

5.3 Final considerations 

To demonstrate the structure of -K^out spectra under study we show in Fig. 3 and Fig. 4 
distributions generated by Herwig [13] for the three Born configurations. Notice that 
the difference between the distributions for each configuration is much more evident 
in the right case since it appears already at DL level. The corresponding numerical 
predictions following from the present theoretical analysis, including the relative 
weight of PT and NP components, will be presented elsewhere [11]. 
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3: Herwig simulation for the total minor (T^ = KJ^^/Q) distribution with 
1 (Durham algorithm [14]). The picture shows the distribution for the three un- 
configurations and their weighted sum. 



We are intended to explore whether the approach we have developed here for the 
analysis of the PT and NP contributions to the i^out distributions and means can 
be extended to more complicated processes and other multi-jet observables. What 
makes one hopeful is the observation that, in spite of technical complications, the final 
results, both in the PT and NP sectors, do possess a simple physical interpretation. 

In particular, we plan to analyze the case with one or two jets in the initial state 
(lepton-hadron and hadron-hadron collisions, respectively). The resulting distribu- 
tion should be sensitive to the structure of the underlying hard matrix element, as 
it is the case for the distributions considered in this paper. Carrying out such a 
programme, and comparing the QCD predictions with experiment, should provide a 
deeper understanding of the hadronization physics. 

Acknowledgements We are grateful to Gavin Salam for helpful discussions and 
suggestions. 
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Figure 4: Herwig simulation for the narrow minor (T^ = K^^/Q) distribution with 
Vcut = 0.1 (Durham algorithm). The picture shows the distribution for the three underlying 
configurations and their weighted sum. 



A. Soft radiation resummation and PT distributions 

The resummation of the PT contributions to the distribution Sj^'^ to SL accuracy 
has been performed in [1]. To achieve such an accuracy one needs at two- loop 
level. The recoil momenta qa need to be taken into account in the kinematics, but 
they can be neglected in M^^. The essential point for the resummation of the soft 
radiation is the factorization of ^ and of dHn^'^ , the soft multi-parton distributions 
and the phase spaces (2.10) and (2.11). In this Appendix we recall the PT relevant 
points needed to analyze the NP corrections in S(ii'out) which emerge as an interplay 
between PT and PT contributions. 

Before turning to Tj^^{Kout) we recall the structure of the soft matrix element 
factorization at two-loop level. We start with the one loop expression, without virtual 
corrections, given by 

n 

Mis ^ l[Ws{h), (A.l) 

i=l 

where Ws is the distribution of soft gluon radiation off the hard three-parton antenna 
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in the configuration S. For instance for 6 — 3 one lias 

Wsik) = — I Wi3 + W23 - Jp^l^ 



a. 



PaPb 

witfi kab,t the invariant transverse momentum of k with respect to the a&-dipole. In 
general we can write 

Ws{k) = J2 C^^^abik) , Cif = Cl? = ^ , C!^^ = ^, c^a,b. (A.3) 

a<b 

This expression can be generalized at two-loop by considering the emission of a pair of 
soft partons from the three-jet configuration [12]. The two-soft emission distribution 
has the same form as (A. 2) and this justifies the exponentiation of the two-loop soft 
emission distribution. As shown in [5], the two loop corrections to the soft matrix 
element enter at PT level simply by properly setting the argument of the running 
coupling in (A. 2) at the dipole invariant transverse momentum. Since for the PT 
analysis at SL accuracy we can neglect the recoil in the soft matrix element, we set 

Pa = Pa in Wab- 

The factorization of the phase space dHn^^ is achieved by using the Fourier and 
Mellin representations for the delta and theta functions in (2.10) or (2.11). Each soft 
parton ki enters with a source u{ki) which, at PT level, must satisfy the condition 
u{ki) — > 1 for ki coUinear to Pa- In the next two sections we recall the form of 
the sources, the radiators and the PT distributions for the and K^^^ cases in 
succession. 

A.l -ft'Jut case: source and PT radiator 



^'^^'^'"""^ -Qe-k..! -Qe-Hfc..l 




27rii/ 

a=l 



n 



u^£he-'-^^^ TT ....... 



e 

R 



(A.4) 



2. n 

u 



e 



n 



27r 

D 



with qf^ = qi^ ^?(±gij/). We recall that R ^ Ci + C4, U ^ Ci + C2 and D^Cs + C^ 
represent the right-, up- and down- regions of the phase space with Q the regions in 
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the four quadrants. Each emitted soft parton contributes with the source (for the 
PT case) 

4 

uP{k) = J2^eik)uf\k), (A.5) 

£=1 

where Q£{k) is the support function for k & and 

uf\k) = e-''{l'=-l-^'=^2fe4 ^ uf\k) = e-''{l*^-l-^^3'=-> , (A.6) 

Pi = bix , 7 = biy , Pl2 = Pi + P2, Pl3 = Pi+ Pa- 

The source uP{k) depends on kx and ky for A; e Ci or C4, and only on k^ for k e C2 
or C3. We obtain (2.12) and (2.13), where 

The source (A.5) satisfies the condition uP{Q) — 1 and therefore constitutes a suit- 
able source for the PT analysis at SL accuracy where recoil can be neglected. 

By using the factorization of soft gluon emission (A.l) and of the phase space 
(2.12) we can resum (2.9) and obtain 



n i n=0 i 

= exp jy [dk] Wsik) uP{k)^ , [dk] 



n=o " i (A.S) 

d^ 

TTUJ 



where, in the r.h.s we removed the integration dfi^, see (2.13). 

Next we include the virtual corrections by subtracting unity to the source and 
obtain the PT radiator for the total distribution 

J [dk] Ws{k) uP{k) ^ - j [dk] Ws{k) [1 - 4°^^)] = -^r(^, Pa, 7) ■ (A.9) 
This expression for the PT radiator of E^^ 

7^^(^,/i^a,7)=E^a?^«'- ^»''= f — Wabik)[l-uPikx)], (A.IO) 

is valid at two- loop accuracy provided one uses the correct argument of ag. One has 

nk! 



^al,{k) = , (A.ll) 



t,ab 
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with the couphng defined in the physical scheme [10] . This radiator has been evalu- 
ated in [1] and is given in (2.23). 

The Mellin moment cr|'^(i/) is given by the d/x^ integral over the three "Sudakov 
form factors" associated to the three hard partons. Since the PT radiator does not 
depend on the recoil, we can freely integrate the function over qiy and Qax and 
get 

4 /. 

+ 



(1 + /5|)(1 + /3|) V(l + /5?2)H7 + e) {l + f3l){t^ + e 

To obtain the Mellin moment we need to integrate over the Fourier variables. Con- 
sider the integral 



TT 



oo ■ 



with Ri = Ri(z/, /3i2, /9i3, 7) given in (2.23). The symmetry of the 7-integrations 
implies that Z is a function of the sum /?23 = (^2 + (^3- The radiator Ri has an 
imaginary part; however, due to the symmetry 

Ri(i^,/3i2,/3i3,7) = Ri(i>,/3i3,/3i2,-7) ={Ri(i^,/3i3,/3i2,7)r , (A.14) 
and to the properties of the /32, /Sa-integrations, the final result is real and reads 

f°° dd 

^(^,/323)= y^ ^(i + ^2) S{v,(5,(5,,), (A.15) 
with S given in (3.6). The real and imaginary parts of Ri are given, respectively, by 

(p, /?', 7) = ifi) + \Ri ifi') , 

r\ (A.16) 

R«(P, P, (5\ 7) = I^M Bi(3, 7) - Kifl') B{P', 7) , 

where 



^^^^(1+^)2 + ^2^ ^'^^^(1 + ^)2 + ^/2^ ^^-^7^ 

and where the function B{P, 7) is given by 

WR ^ r dx , (l + x7)2 + /32 



+ 7)' + 



In a; 



(A.18) 



/o 27r(l-a;2) (a; + ^)2 + 2;2^2 
This function has the following limits 

5(/3,7) ^ Ofor l,7«/9, 5(/?,7) ^ for l,/3«7- (A.19) 
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A. 2 -f^^ut case: source and PT radiator 

Wc factorizc the phase space by using MeUin and Fourier representation and we 
obtain (2.12) and (2.13) with 

V^i(3, qi^) = e-^(l«i=^l-^«i=^^) , (A.20) 

and the source for the soft parton emission 

uf{k^) = ii(°)(A;) QRik) + QLik) , u^°\k) = e'"^^''^^ - "^''^^ , (A.21) 

where QR/L{k) is the support function for k & R/L. This source satisfies the condi- 
tion u^^^ (0) = 1 and therefore is a suitable source for the PT analysis at SL accuracy. 

Proceeding as before, to SL accuracy, the PT radiator for the distribution 
is given by 

= E ' = / ^«''(^) [1 - , (A.22) 

Wahik) given in (A. 11). To SL accuracy the i^^^'^^^iator is given by (2.40) with 
computed in [1]^. 



B. NP correction to the radiators 

In this appendix we compute the NP correction to the radiators for the K^l^ dis- 
tributions. Wc follow the method used in [5] to compute NP corrections to various 
coUinear and infrared safe observablcs in two-jet events. 

We start by recalling the relevant points in the calculation of the NP corrections 
to the radiators for the two-jet CIS observables which we need to generalize to the 
present analysis. We focus our attention on the broadening observable [7] which 
takes contributions from the transverse momentum of the emitted particles. 

B.l Recollection of two- jet case 

Kinematics. One introduces two opposite Sudakov momenta P and P, with P 
along the thrust axis. For a (secondary) parton k of mass m the Sudakov variables 
are given by 

P = |(1,0,0,1), P = |(l,0,0,-1), k^aP+(3P + kt, a(3='^^, (B.l) 

with kt orthogonal to P and P. The right region R is defined as 

keR: a>/3 =^ a> + . (B.2) 

^This scale depends on the configuration index 5. In this reference we denote this scale for the 
configuration 6 as Q^^ 
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To compute NP corrections to the broadening one needs to take into account the 
recoil of the primary quark and antiquark momenta for which we have a similar 
Sudakov decomposition. 

Splitting of the radiator in three pieces. The full radiator is a CIS quantity. 
One rearranges [5] the two-loop radiator as a sum of three CIS terms which have been 
called naive, inclusive and non-inclusive pieces. The naive piece is the contribution 
for the emission of a single massive gluon together with the corresponding virtual 
corrections. Its form is given by 

where a is the Sudakov variable, k the transverse momentum of the soft massive 
gluon with respect to the hard emitting parton (including the hard parton recoil). 
The integration is over the right hemisphere (B.2) for the massive gluon. The effective 
coupling aes{m'^) is given in terms of the physical coupling in [8] through a dispersive 
representation. 

Massive source. As already stated, the source u{k) takes into account the recoil of 
the two hard partons and is then given in term of the relative transverse momentum 
K of the massive gluon. One is free to chose any prescription to implement the mass 
dependence in the source, provided that u{k) reduces to the massless case for m ^ 0. 
Following [5], we assume that the massive gluon source is given by the expression of 
the massless gluon source in which we make the replacement 

K, VkP' + m? . (B.4) 

Milan factor and universality. Using this prescription for the massive source, 
one has that the naive piece is the only one that needs to be computed since the 
remaining pieces (the inclusive and the non-inclusive pieces) contribute in a universal 
way simply by rescaling the naive piece by a universal Milan factor. 

Linearization of the sources. The leading NP correction coming from (B.3) is 
obtained by taking the linear contribution in m and k of the integrand, that is the first 
term in the expansion of [1— ■u(A;)]. This fact has important implications on which are 
the relevant kincmatical and Fourier variables affecting the leading NP corrections. 
They arc the variables which enter the expression of the linearization of [1 — 'u(/c)] 
and do not vanish upon the integration over the azimuthal angle of R in (B.3). In 
the case of broadening [7] the Fourier variable 6, conjugated to the momentum k 
does not enter the leading NP correction to the radiator. Its hnear contribution to 

— * 

[1 — ii(A;)] is proportional ioh ■ k which vanishes upon azimuthal integration. 
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B.2 Sudeikov variables and soft distribution for the a6-dipole 

We now return to our three-jet case and we start to discuss the kinematical variables. 
The radiator is given by the sum of three dipole contributions Tab (for the PT case 
see for instance (A. 10)). To treat each dipole contribution Tab by the same method 
used in the two-jet case for (B.3) one needs to perform a Lorentz transformation. 
Consider the e+e^ laboratory system defined by the Sudakov variables introduced in 
(B.l) with P along the thrust axis. In this system the two massless partons Pa,Pb 
initiating the a6-dipole emission are not in opposite directions. In order to follow 
the method of [5], we use the fact that the distribution Wab is Lorentz invariant and 
we analyze rab in the afe-dipole center of mass (cm.). There are three dipoles and 
then there are three different cm. systems. Given a dipole ab, we introduce the 
a6-dipole cm. and express the distribution Wab{k) in these variables. We then apply 
the two-jet case analysis to our case. We describe the new variables. 

Dipole ah center of mass frame. Consider the laboratory frame in which the 
Born momenta Pa, Pb of the a6-dipole are given in (2.3). We go from this frame to 
the a6-dipole cm. frame by a rotation in the y2;-plane, a z-boost and again a rotation 
in the yz-plane which takes and P^ along the 2;-axis. We denote by P^, Pl,p'^,p'^, 
and k'^ the parton momenta in this frame with Sudakov components 

p: = ^(i,o,o,i), pi 

p'a = KbK + KbPb + , 
Pb = <aPa + B'baPl + Q'bt , 

The transverse components are orthogonal to P^ and P^'. All Sudakov variables are 
soft except for A'^^f^ and P^,^ which are close to one. The x-components are not afi^ected 
by the Lorentz transformation which is in the yz-p\a.ne, so that the components 
entering our observable Xout are the same in both frames 

We denote by Fab the forward region of secondary parton jj^i in the a6-dipole cm. 

KeFab-. a\>l3[ ^ a[>^. (B.7) 

The Sudakov variables (B.5) are different for different dipoles. We do not need to 
add an index to specify the dipole since these variables will be used only in the 
calculation of the specific a6-radiator. 



%(1, 0,0,-1), Ql,^2PaPb, 



A' ~ 1 

^ab — ^ 

R' ~ 1 

^ba — ^ 



it 



a'-ff- - 



ab 



(B.5) 
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Soft distribution in the a6-dipole cm. The (Lorentz invariant) a6-dipole soft 
distribution Wab is given in terms of kab,t, the "relative transverse momentum in the 
a6-dipole frame". By using the Sudakov variables (B.5) in the a6-dipole cm. frame 
we have 



k1 



k 



(B.8) 



where we have neglected the large components A'^^ ~ B^^ ~ 1. 
We split the dipole distribution Wab into two pieces: 

Wab = <6 + , (B.9) 

where has support only in the forward region Fab, 

wl,{k) ^ wab{k) ■ e(a'-(3') ^ ^ ^ ^ _ ^.^-..^ _ (B IO) 

The "backward" piece w^j^ is similarly defined. Notice that Wab{k) is invariant, but 
the splitting in the two regions is performed in the a6-dipole cm. 

In three-jet events the soft distribution W^^\k) is given by a sum of dipole soft 
distributions Wab{k), see (A. 3). We decompose Wab into the two contributions w^b 
and w^fj as given in (B.9) and we have 

iy(^)(A;) = 5^Ci?<,(A;), 

(ab) 

where the sum is extended to the six ordered pair {ab), i.e. ab = 12, 13, 21, 23, 31, 32. 
Each term contributes only in the region Fab (recall that has a support only in 
the forward region (B.7)). 

B.3 Source with recoil for the K^y^^ distribution 

As discussed above, the leading NP corrections to the radiator are obtained from 
the first order expansion of [1 — liT(^)] in the momentum components. This implies 
that, as in the broadening case, the leading NP radiator does not depend on the 
four Fourier variables and 7. Their first order contributions in the expansion 
of [1 —UT(k)] are linear in the momentum components and then vanish upon the 
symmetric phase space and Fourier variable integration. We have that only the 
Mellin variable enters the leading NP radiator in the form ulki^l. The conclusion is 
that in the phase space dHn^^ in (2.10) and (2.11) we should take into account the 
NP gluon only in the theta function giving the upper bound i^^out • On the contrary, 
the NP soft gluon does not contribute to the determination of the event plane. The 
delta functions which fix the event plane should be taken into account only in the 
PT contribution. 
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We now analyze how one can take into account in the sources the recoil coming 
from the observable K^y^^. As in (A. 8) we start from the sum 

T.l^J^HnU^dh]Mls = J^^^^^ Jd^^expU[dk]Ws{k)uP{k)\. (B.12) 

Using the expression in (B.ll) for the soft gluon distribution we can write 

exp I [ [dk]Ws{k)uPik)\ = Hexp {^1? [ [dk] u^ik)] . (B.13) 

^ {ab) ^ -J J 

From (B.IO) one has that w'^ij{k) is coUinear singular for 

K^k[-a'q'^,^Q, k,-a'qa,^0, (B.14) 

(with the prime variables defined in (B.5), the afe-dipolc cm.). From the form (A. 5) 
of these source we have u^ik) 7^ 1 in the limit (B.14). We then need to construct 
a new source uxik) which, in this limit, satisfies the condition [l—UT{k)] — > 0. This 
new source must include the hard parton recoils. Actually we need to consider the 
recoil only for the contribution to the source associated to the theta function fixing 
the observable -fi'^t- To do this we start by expanding the functional in (B.13) 

exp I J [dk] ws{k) uPik) I = n E n ^^'^V [^^^^ """^^^^^ ""^^^^^'^ ■ 

^ (ab) riab " ' «=1 

Consider a given term in this expansion. For any given ordered pair (ab) there are 
Uab soft gluons which are associated to the distribution w"^ and are then emitted 
in the region Fab- For this term in the expansion (B.15), the corresponding theta 
function fixing the observable Kj^^ gives 

Kout = E + E l^'^l ^ E + E 1 E l^^^l f • (B.16) 

a i a (ab) KieFat, J 

Recall that the x-components are the same in both the laboratory and the a6-dipole 
cm. frames. We then introduce the subtractions 



ix \-<^'^Qax\) ? , 
a (ab) KieFab ) 

qax = ( 1 + E E ^i] 

\ b^a ieFab / 



(B.17) 



(a^ is defined in the afe-dipole cm.). Here we have neglected corrections of second 
order in the soft parameters. 
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Now we reconstruct the expansion (B.15) by using (B.17) instead of (B.16). 
While (B.16) gives (B.12), using (B.17) we obtain 



n i \ a=l 



Hexp Ci? / [dk]w2,{k)u,,{k)\ , 

{ab) ^ 



(B.18) 



where Uab{k) is a source associated to the ordered pair (ab) given by 

Uab{k) = ufik) ■ e^"'l«»-l = e-''(l*^-l-"'l«»-l) • UT{k) , (B.19) 

where urik) is the part of the source depending on the Fourier variables /3a, 7 (see 
(A.6)). The factor g-'^d'^-l-"'!^'^-!) , which is the only one relevant for the NP correc- 
tions, is equal to one at the singular point (B.14). We have then that (B.19) is the 
form of the source we are searching for. 

The Jacobian for changing integration variables from to qax is given by 

a \ b^a i&Fab / {ab} i^F^b 

This can be taken into account as a correction to the soft distribution w'^i^{k) — > 
{l — a')w'^f^{k). This correction corresponds to including the hard part of the splitting 
function, see [5]. Since in (B.18) the rcscalcd momentum qax is just an integration 
variable, in the following we shall rename it qax- 

B.4 NP radiator for the K^y^^ distribution 

The radiator for the K^^-distribution, given in terms of the dipole components (B.9), 
reads 

^^ = E^a?^«^' = r:,^ f^w:,{k)[l-ua,{k)], (B.21) 

T^b J 

where the integration region is restricted to the Fab region. 

The distribution for the emission of two soft partons in the three-jet event (see 
[12] and [1]) can be represented as the sum of the two soft parton emission from the 
three dipolcs. This fact can be used to obtain the two-loop PT radiator as sum of 
the two-loop radiator for each dipole. At this point we can use all the results of the 
two-jet analysis. At PT level this gives rise simply to the identification of the running 
coupling scale for each dipole radiator, the invariant dipole transverse momentum, 
and the scheme. For the analysis of the NP corrections we follow the two-jet case 
method recalled above. Introducing the coupling according to the dispersive method 
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[8] we write the two-loop dipole radiator as a sum of three pieces, the naive, the 
inclusive and the non-inclusive pieces. The naive piece of the a6-dipole is given, as 
in (B.3), by 



r 



/, 9 / 9n —d f (Pk f da' , ,,,, 

dm'a,sim')-— / — [1 - Uabik)] , B.22 

Qab 

with Q;efr(m^) the effective coupling discussed in [8]. The integration region is re- 
stricted to the forward region Fab defined by a' > k[/Qab- This has been replaced 
by a' > y/ + m'^/Qab, since, at the lower limit of a' (large angle emission) the 
correction is of second order in the NP-parameter (see [7]). 

As discussed in the Appendix B.l, the leading NP correction coming from (B.22) 
is obtained as follows: 

• in the integrand we take the linear contribution in m ~ Aqcd <S \qax\- 
Therefore, for the massless gluon source (B.19), we take the first term in the 
expansion 

[1 - Uah{k)] = u{\kx\ - a'\qax\) + ■ ■ ■ k^c = K cos (p + a' Qax , (B.23) 

where (p is the azimuthal angle of k. The neglected terms contribute to the 
higher order power corrections; 

• for a gluon with mass m we make the replacement (B.4). This ensures a 
simple treatment of the other (inclusive and non-inclusive) pieces of the two- 
loop radiator; 

• we take the NP piece 5aef{{m'^) of the effective coupling. 
The leading NP correction to the naive piece is then given by 



Qabi^''^ + m'^) ^ J ^ J^r^r:^ (^Wk^ + iti^ cos (p + a' qax\ - a'\qax\) 



Qab 

To evaluate this integral we introduce the rescaled variable 

Oi'\qax\ 
V — . = 

V k'^ + w? 

and find 

dv 



(B.24) 



Q„5(«' + m^) = i/V^c^ + m2 r ^ ! —{\cos4>-v\-v). (B.25) 

J-TT J\Qax\ V 

Qab 
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Notice that at the upper hmit a' = 1 one has v = which tends to infinity in 

the hmit we are considering. On the other hand the contribution with f > 1 vanishes 
due to symmetry of ^-integration. Setting q^^ — > at the lower hmit, we find 

^ab{i^'^ + m^) = ck^^^ \/k^ + w? In , ck^^^ = - , (B.26) 

with ck^^^ the characteristic number for the Xout observable and which takes into 
account the small- angle contribution, is given by 



lnC = I / " ^ I ' ^{\cos<P-v\-v -\cos<P\] =lii2- 2. (B.27) 



We get 



KbiQax) = In j—^ , A = ck^^^ / — m . (B.28) 



m- 



For the "inclusive" and "non-inclusive" pieces the analysis is the standard one and 
we have that the full NP correction is given by (B.28) by rescaling the NP-parameter 
by the Milan factor }A. 

The NP parameter A'^^ has dimension of a mass, and this imphes that the NP 
correction is proportional to 1/Q. It can be expressed in terms of the integral of the 
running coupling over the infrared region 

Mf^i) = / . (B.29) 

Jo 

This parameter aol/^/) is the same entering the two-jet shape variables 1 — T, C, 
B and jQ^. After merging PT and NP contributions to the observable in a 
renormalon free manner, one has that the distribution is independent of /// and one 
obtains 

A^^ = ck^^^ ^^i^i |«o(/^7) - - /3o^ (^^ ^ ^ ^ ^ } ' ^ "ms(<3) ■ 

(B.30) 

The term proportional to K accounts for mismatch between the ms and the physical 
scheme [10] and is given by 

K = Ca\ \--n(, /3o = . (B.31) 

It may be useful to connect this parameter entering the shift of the i^out distributions 
to the analogous one for the thrust distribution. Defining r = 1 — T, we have 

da , , da^^ , ^ ^ * ^ Cr A^^ ,^ 

— T = — — r-A, , K = Cp ^ , B.32 

dr dT ck^^^ 
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where Cp enters due to the fact that the two-jet system is made of a quark-antiquark 
pair, Ct = 2 and ck^^^ =2/11 are the characteristic numbers associated to the phase 
space integration of the relative observables. 

We can now reconstruct the NP correction to the full radiator in (B.21) and we 

find 

NP i„ CQ12 , ^ CQ12 , ^ CQ23Q13 

Qi2\Q3x\ 

iNP/^ i„ ''Via , ^ i„ '.Vis , ^ , CQ23Q12 



Cf 


In 


CQ12 


+ Cf 


In 


CQ12 


+ Ca 


In 






\(llx\ 






\q2x\ 






Cf 


In 


CQ13 


+ Cf 


In 


CQ13 


+ Ca 


In 






\Q1x\ 






\Q3x\ 






Cf 


In 


CQ23 


+ Cf 


In 


CQ23 


+ Ca 


In 






\Q3x\ 






\(l2x\ 







5nM = 

5'R2{q) = vX^"" ( Cf In + Cp In + Ca In l^ipii ) , (B.33) 



NP/^ 1„ ''V23 , ^ i„ ^V23 , ^ , CQ12Q13 



Q23|?la 



We find then that the NP scales are given by the hard scales Qa in (2.22) rescaled 
by a factor = C (see (2.21)). 

B.5 NP corrections to the -ftT^ut distribution 

We proceed as in the previous case. Here the simplification is that we should be 
careful only with the collinear singularity with respect to pi which is contained only 
in the distributions w\^{k). 

First we consider the modification of the sources. We must include the recoil 
only in the source un,{k) associated to the distribution ^^^(/c). We have 

uib{k) = QLik) + QRik) e-^(l^-l-"'l«i-l) e^^'^- . (B.34) 

For the other sources Uab{k), associated to the distribution w'^j^{k) with a 7^ 1, we do 
not need modifications: Uab{k) — u^^\k). The modification in (B.34) is accompanied 
by a change of variable from qi^ to qix obtained by a rescaling factor similar to 
(B.17). However since qi^ is an integration variable we shall rename it qi^- 
Then we take the linear contribution for the massive gluon source 



[l-uib{k)] =^ u (^\V li^ + w? cos(/) + a'qix\ - ©r(/c) , 

[1 — Uab{k)] =^ + m?\ COS0I 6i^(A;) , a^l. 



(B.35) 



Since the distribution has a singularity for the soft gluon collinear to pi, the NP 
radiator has a logarithmic behaviour in qi^- We have 

(QNP 

SUsiq) = i/A^^Cf ) In . (B.36) 

The In \ qix\ part is contained in and rjg since the corresponding distributions are 
collinear singular in the integration region. The remaining pieces contribute only to 
set the NP scale. To obtain the NP hard scale Q^^ we need to analyze the individual 
pieces. 
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To see the singular behaviour in (B.36) consider the NP correction to rj^. At 
two loop, the naive term of is given in terms of the effective coupling as follows 

t\, = [ dm'a^^im')-^^ [ f u,,{k)] . (B.37) 

Qlb 

Taking the linear part of the source and NP part of the effective coupling, we have 

/—d f dK^ 
dw?5ae^{w?)—^ / ^— 5- fiib(fi;^ + m?) , 
drw^ J K + 

/TT ^0 j'l (B.38) 
— / ( I COS0 — v| — w ) 9r(/c) , 

27r Jkaxi V 

Qla 

where we need to express the condition in 6r in terms of the variables in the 16-dipole 
cm. system (see (B.5)). 

B.6 Recoil integration over the NP radiator 

For the K^^^ case we have 

" (B.39) 
^ 7) 1 1 - ^A^'' E ^« [ l^(^Qr) + Xa(/?, 7) ] I , 

with V'^ and /(/5,7) given in (A. 7) and (A. 12) respectively, X2 = x(/32), Xz = xiPs) 
(see (3.2)) and 

(! + /?!) (! + /?!) Ui + /??2)H7 + e)"'(i + /??3)(^7 + e)y'- ^ • 
For the K^^^ case we have 



where is given in (A. 20). Both in (B.39) and (B.41), corrections are of order 

(^NP)2_ 



C. The NP functions fa 

C.l /i 

To compute /i we first split the function x(/3) in the two pieces in (3.2) and evaluate 
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the two contributions separately. We have /i(z^) — /{(z^) + fi{^) where 



(C.l) 



°° a=2,3 



The evaluation /{ is straightforward since all three /3-integrals are rapidly convergent. 
For large v we can factorize the three radiators and the remaining piece is a function 
of the SL variable r'{v, Q) 

/(H = ne-«"^'^-C;(z.), CU^)=i?;m(K(^)}) . (C.2) 

a=l 

For i/ = 0, we set to zero the radiator and, due to the integration property of x' in 
(3.2), we have /{(O) = 0. 

To evaluate /f we observe that the (32- and /33-integrals are rapidly convergent, 
so that the characteristic values of these two variables arc (32-, = On the 

contrary the /9-intcgraI is logarithmic so that the characteristic values of /? are large. 
From (A. 16), (A. 17) and (A. 19), in the region 

P2, Ps<.f3c^f3T, (C.3) 

we have 

Rf)(P,/3,/3T,0) ~ R,(P(3) Rt\p,(3,(3T,l) ^0, (C.4) 

(see (A. 16) and (A. 19)). Therefore, in S one can factorize the ^1 radiator and the 
remaining factor is a SL function 

S(u, (3, (323) = e-^^^") • { 1 + R[(P) B, ((3, (323, K(m , P = , (C.5) 

with Bi vanishing in region (C.3). The first contribution gives the leading contribu- 
tion to /i. We can factorize the P2- and /Ss-integrals and, changing variable from /3 
to p we obtain 

(C.6) 

where we have factorized the three radiators, the functions J^a are the single jet 
functions defined in (2.44) and the function Ei is introduced in (3.9) and discussed 
in detail in Appendix D.l. 
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Consider the second term contribution in (C.5) to It will be denoted by /f . 
Since Bi vanishes for large f3, all three /3-integrals are rapidly convergent. As for f[, 
we can factorize the three radiators and the remaining integrals give a SL function 

/T(^) =fle-^^^'^-C';{u) , C'liu) = R[{u)A[ {{K{iy)}) . (C.7) 

a=l 

Here again we have that C" vanishes for R[ — > and i/ — > 0. Assembhng all pieces 
(C.2), (C.6), and (C.7), the function /i(i^) is given by 

3 

/iW = n^""""^"^ {M^) J'siu) E^{P) + CiH} . (C.8) 

a=l 

Here Ci(i/) = C'llu) + C'liy) is a SL function, which vanishes for E!^ —>■ and for 
t/ -> 0. . 

To conclude we compute the first loop contribution proportional to R'^ to S. We 
have 



y 1 + /?r TT Jo 7 

(C.9) 

The first two terms are the 7 = contribution. The last term has been evaluated in 
[1] and one finds 

— ImRi(t/,/3,/3T,7) = — — -In^ (C.IO) 

7 2 y^l + /J^2 



2 r 



Thus the one loop contribution is given by twice the second term. It vanishes for 
C.2 A 

As before, to evaluate /2 we first split the function x{l32) in the two pieces in (3.2) 
and compute the two contributions. We have /2(z^) = f2i^) + /2 ('^)) where 

^^^^^ ^ Xoo ^3 .H + P^J L + ^ ^^'^ ' 



^'^'^^"70 1 + loo vr(l + /3|) I 



7r(l + /32) ■ 

(C.ll) 



The analysis of /2 is similar to the one of f[ performed before and we get 

3 

m = ne""""^'^ •^2(^^) , c'^ii^) = i?^(p)^2({<(p)}) , (C.12) 



a=l 
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where C^i^u) is a SL function vanishing for R2 ^ and u ^ 0. 

To evaluate /2 we observe that the ^3- and /3-integrals are rapidly convergent, 
so that the characteristic values of these two variables are ^,^3 = O {1). On the 
contrary the /32-integral is logarithmic so that the characteristic values of P2 are 
large. Therefore, to evaluate /2 we need to obtain S{p, /3, /323) in the region 

/?3 « . (C.13) 

We show that here one has 

S{u,P,P23)=eM'^) .{l + R[{u)B,{u,P,P,s)}, (C.14) 

where B2 is a SL function of u and rapidly vanishes in the region (C.13). 

To obtain this result we observe that in this region we must consider also the 
second term of S in (3.6) with the 7-integral and the imaginary part of Ri. This 
integral is logarithmic, as the /32-integral, so that the characteristic values of 7 are 
large compared to P,f33. On the other hand, from (A. 16), (A. 17) and (A. 19), we 
have R^*^!/, /3, /?t, 7) — >^ for (3t <C 7. We conclude that the region 



P, Ps « 7 « /3t ^ /92 ^ P23 , (C.15) 

gives the leading logarithmic approximation of f^- We start from this approximation 
which we shall later improve to account for the first subleading correction O {R[). 
In this region, from (A. 16), we have 



R«(z7,/5,/5t,7) ^ ^i?i(^7) + hR^m) , 

71 



and, for the two terms of S in (3.6), we have 



(C.16) 



g-Ri(P,/3,/3T,0) _ g-5iii(P/32)g-5«i(i'0+^) ^ 

rtLime-^^i^'^'f'-'^^ ~ e-h^^^'^^) r"^e-^«^(^^) sinf-i?;(i.7)V ^^'^^^ 
Jo 7 A™ 7 V4 V 

where jm = 0(^2) and 7^ = O {P). To evaluate the contribution from the imaginary 
part we first consider 



sm 



(^^'1(^7)) = + (C.18) 



The terms in the dots can be neglected since they contribute to the integral with 
terms of order ccg. To show this consider the contribution to the 7-integral from the 
first term which is of order (cts In 7)^ for large 7. Since the 7 integration is logarithmic 
and the radiator is of order Q;sln^7, the result of the integral is of order ctg. 
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Taking the first term in (C.18), we find the leading logarithmic contribution 

Here we used jm ~ /32- 

In order to keep under control the subleading O {r') correction we need to fix 
precisely the lower integration hmit 7^ in (C.19). To do this it suffices to calculate 
the exact 7-integral in the first order in r' given in (C.IO). Combining (C.19) and 
(C.IO) we have the result 

TT 7o 7 

(C.20) 

Summing this with the real part contribution in (C.17) we have the cancellation of 
the leading term e-5^i(^*)e-5«i(^v/^) and we obtain 



5(L/,/3,/323)^e-^^(^), p=P^l + /3|. (C.21) 

This result shows that if one factorizes in S the radiator e"^^^^\ the remaining piece 
is a SL function and one obtains the final result in (3.15). 

The leading piece to /2 is obtained from (C.21) and is given by 



3 



(C.22) 

with E2 given in (3.17) and given by (2.44). 

The second term in (C.14) (with B2) gives, after factorizing the radiators, a SL 
function since all three /^-integrals are rapidly convergent. One gets 

~fi^'{v) = fl e-^^^'^ ■ CM , Ci{v) = R[{u) A',{u) , (C.23) 

a=l 

with A2 a SL function. 

Putting together the various pieces, (C.12), (C.22) and (C.23), we get the final 
result 

3 

/2M = n^'"""^'^ • {^3{^)E2{i^)+C2{i^)} , (C.24) 

a=l 

where €2(1^) — R'2{i>)A2{u) + R'i{p)A'2{u) is a SL function which vanishes for r' ^ 
and i/ — > 0. 
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D. The E functions 

For completeness we recall here the results of an Appendix of Ref. [7] and extend 
them to the present case. 

D.l El 

We consider first the function £'i(i>). For the function Ei[i') we have the expansion 



J u P 



with R[,R'(,Ri the first, second and third logarithmic derivatives of Ri evaluated 
at P and 

N{t) ^ A Hdx e--^+*^ , X{t) ^ = r dzz e— 2*v^ , (D.2) 

V^'' Jt ° Jo 

with the following behaviours 

(D.3) 



In the region t — R[/y'2R'[ ~ ^/a^\n{lyQ) 1, we obtain 



On the contrary, for t — R'i/^/2R'( ~ ^/^i^ln{uQ) -C 1, we find 

2^/Cias{Q) D<^i 

Notice that the term O implicitly includes all corrections O (y^) : O (y^ln z/) 

and O (i/a^lnP^). 

D.2 Ea for a = 2, 3 

For the function Ea{v) we obtain the following result: 

f) 



50 



where again i?^, i?^', are the first three derivatives of Ra evaluated at v and the 
functions N and X are defined in (D.2). 

As in the previous case, in the region [R'^ + i?^)/ iy2(^R^~+~R^ ^ 1, using the 
expansions in (D.3), we get 

^a(^) = o.... ^ o,... + ( , i ] . (D.7) 



On the contrary, for (i?; + R'^) / ^2{R!{ + R'^) < 1, we obtain, up to contributions 



^a(P) ^ ^ ---- C,ln(pgr)+C„ln(pgr^) + ^ ■ (D.8) 
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